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ABSTRACT 


The  coupling  was  made  between  a  structural  analysis  code 
(VEC/DYNA3D)  and  an  underwater  shock  analysis  code  (USA) .  The 
coupled  computer  program  (USA/DYNA3D)  was  verified  using  a 
set  of  benchmark  problems  which  had  known  analytical 
solutions.  The  benchmark  problems  were  elastic  analyses  of  a 
spherical  and  an  infinite  cylinder  subjected  to  a  plane  wave. 
The  comparison  between  the  numerical  and  analytical  solutions 
was  very  good. 

An  underwater  explosion  test  was  performed  with  an 
aluminum  cylinder  subjected  to  a  far  field,  side  on 
explosion.  A  pre-shot  calculation  using  USA/DYNA3D  determined 
critical  locations  to  measure  both  axial  and  hoop  strains  of 
the  cylinder.  After  the  experiments,  a  post-shot  calculation 
was  undertaken  using  the  free  field  pressure  obtained  from 
the  physical  experiment.  The  numerical  results  obtained  using 
the  elastoplastic  analysis  were  very  comparable  to  the 
experimental  data  except  for  some  positions. 

A  series  of  numerical  experiments  were  performed  to 
determine  the  cause  of  the  difference  between  the  numerical 
and  experimental  results  at  some  positions.  It  was  found  from 
the  experimental  data  that  there  might  be  some  rotation  in 
the  cylinder  relative  to  the  explosive.  Considering  the 
rotation  in  the  numerical  model  improved  the  comparison.  In 
addition,  numerical  sensitivity  analyses  were  undertaken  to 
determine  the  importance  of  various  physical  and  numerical 
modelling  factors. 

Finally,  this  study  showed  there  were  three  major 
response  modes  of  the  cylinder  subjected  to  a  side  on 
explosion.  They  were  the  accordion  mode,  breathing  mode  as 
well  as  the  whipping  mode.  Large  plastic  strains  occurred  at 
the  center  of  the  cylinder  on  the  reverse  side  to  the 
explosive  and  near  the  ends  of  the  cylinder*  on  the  near  side 
to  the  explosive.  The  large  plastic  strains  seemed  to  be 
related  to  the  whipping  motion. 
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I .  INTRODUCTION 

A  research  program  is  underway  at  the  Naval  Post  Graduate 
School  to  study  numerical  modeling  of  ship  structures 
subjected  to  both  near  and  far  field  underwater  explosions. 
This  program  is  expected  to  improve  the  understanding  of 
factors  affecting  the  reliability  of  numerical  models.  In 
addition,  it  will  provide  insight  into  the  dynamic  response  of 
surface  ship  and  submarine  hulls  and  the  physics  that  lead  to 
failure  when  a  hull  is  subjected  to  an  underwater  shock  wave. 
The  current  study  centers  around  simple  cylinders  constructed 
of  a  homogenous  material.  Future  studies  will  include  more 
complex  materials  and  structures  as  experience  increases  and 
the  reliability  of  the  numerical  models  is  proven. 

This  paper  describes  progress  of  the  research  program  to 
date  and  the  expected  line  of  future  research.  Previous 
results  of  this  research  program  were  provided  in  references 
1-2.  This  report  v;ill  describe  the  results  of  two  numerical 
verification  tests  which  were  performed  to  prove  the  validity 
of  a  computer  code  software  link  which  was  developed  for  this 
project.  It  will  also  compare  the  numerical  results  with 
experimental  results  obtained  from  a  underwater  explosion  test 
of  an  aluminum  cylinder  subjected  to  a  side  on  attack.  In 
addition,  analyses  were  performed  to  determine  the  sensitivity 
of  the  results  to  mesh  refinement,  boundary  effects,  rotation 
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from  expected  configuration,  and  use  of  different  types  of 
shell  elements.  Results  show  that  the  computer  code  closely 
models  simple  knovm  analytic  results,  and  can  provide  close 
correlation  to  experimental  results.  Most  of  the 
inconsistencies  between  experimental  and  numerical  results  are 
most  likely  caused  by  uncertainties  associated  with  physical 
model  fabrication  and  the  underwater  explosion  test  rather 
than  a  failure  of  the  numerical  method  to  provide  correct 
answers.  Recommendations  will  be  provided  to  improve  control 
of  future  underwater  explosion  tests.  Finally,  preparations 
for  future  testing  will  be  described  and  recommendations  for 
additional  study  are  provided. 
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II.  NUMERICAL  CODE  DESCRIPTION  AND  VALIDATION 


A.  GENERAL  DESCRIPTION  OF  THE  METHOD 

The  primary  purpose  of  this  study  was  to  match  numerical 
results  obtained  from  a  computer  program  with  experimental 
results  from  an  underwater  explosion  test.  For  this  initial 
study  of  a  side  on  attack,  a  relatively  simple  model  was  used. 
The  model  consisted  of  a  unstiffened,  right  circular  cylinder 
constructed  of  a  homogeneous  material  submerged  in  water.  The 
cylinder  was  modelled  using  a  dynamic  finite  element  m.ethod 
(FEM)  code  called  VEC/DYN.^3D  and  the  water  was  modelled  using 
a  boundary  element  method  (BEM)  code  called  USA  (Underwater 
Shock  Analyser) .  The  linkage  between  the  two  codes  was 
developed  in  1991  at  the  request  of  the  Naval  Post  Graduate 
school  under  funding  provided  by  the  Defense  Nuclear  Agency 
(DNA)  . 

B.  VEC/DYNA3D  FINITE  ELEMENT  METHOD  CODE 

VEC/DYNA3D  [Ref.  3]  is  an  explicit  finite  element  code. 
It  has  been  used  successfully  for  various  types  of  nonlinear 
engineering  problems  since  its  conception  in  1976.  VEC''DYNA3D 
was  selected  for  this  study  for  several  reasons.  First,  as 
stated  above,  VEC/DYMA3D  is  an  explicit  code.  This  attribute 
has  two  distinct  advantages  and  two  disadvantages.  The 
advantages  are  its  relatively  high  speed  and  its  ability  to  be 
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implemented  on  a  relatively  small  stand  alone  engineering  work 
station.  Initial  work  for  this  study  is  being  performed  on 
IBM  RISC  6000  work  stations.  Once  the  USA/DYNA3D  interface  is 
proven  to  be  reliable  and  accurate  and  experience  has  been 
gained  in  the  use  of  the  software,  work  will  begin  on  more 
complex  models  using  main  frame  type  computers.  Therefore  it 
was  important  to  obtain  a  code  that  was  able  to  work 
significant  problems  on  a  small  work  station  and  yet  be 
compatible  with  the  main  frames  expected  to  be  used  in  the 
future.  DYNA3D  is  compatible  with  a  full  range  of  engineering 
work  stations  and  has  been  implemented  on  the  Los  Alamos  CRAY 
computer.  Problems  including  up  to  20000  solid  elements  have 
been  run  on  v;ork  stations  with  16  megabytes  of  random  access 
memory . 

The  first  disadvantage  associated  with  the  explicit 
numerical  code  is  that  the  code  is  not  inherently  stable. 
This  means  that  any  problem.s  dealing  with  time  integration, 
including  the  underwater  shock  probleiCis  included  in  this  study 
must  be  treated  with  care.  Integration  time  steps  must  be 
matched  closely  with  the  size  of  the  elements  in  the  problem. 
This  is  performed  automatically  by  DYNA3D  in  the  stand  alone 
mode.  However,  when  coupled  with  the  USA  code,  this 
automation  is  no  longer  functional.  Incorrect  selection  of 
integrati'^n  time  steps  can  lead  to  significant  oscillations 
and  inaccuracies  in  the  final  solution.  The  second  problem 
associated  v;ith  the  explicit  codes  is  the  mesh  reflection 
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effect.  Non-uniform  meshes  result  in  inaccurate  solutions  due 
to  mesh  reflection.  Two  factors  appear  to  be  important  in 
ensuring  that  correct  solution  was  obtained.  The  first  is 
mesh  size  and  the  second  is  total  mass  of  neighboring 
elements.  Sensitivity  analyses  indicate  that  error  in  the 
final  solution  is  relatively  small  if  neighboring  elements  are 
kept  within  ten  percent  of  each  other  in  size.  This  was  used 
as  a  rule  of  thumb  in  performance  of  this  study.  This  lead  to 
some  inefficiency  in  obtaining  solutions  since  often 
refinement  had  be  performed  over  a  larger  area  of  the  mesh  to 
obtain  a  mesh  independent  solution  than  might  normally  be 
required  in  an  implicit  code.  The  additional  area  means  more 
total  elements  and  a  subsequent  increase  in  computation  time 
to  obtain  the  problem  solution.  These  disadvantages  can  be 
overcome  through  careful  planning.  In  general,  they  did  not 
significantly  overshadow  the  benefits  associated  with  using  an 
explicit  code. 

The  second  reason  for  selecting  VEC/DYN.^3D  was  its  wide 
range  of  available  material  models  and  equations  of  state 
including  the  ability  to  model  strain  rate  sensitivity, 
explosive  materials  and  acoustic  media.  In  addition, 
VEC/DYNA3D  has  a  large  degree  of  interactivity  when  used  with 
the  INGRID  pre-processor  [Ref.  4]  and  TAURUS  post-processor 
[Ref.  5] .  Changes  can  be  entered  with  relative  ease  using  the 
pre-processor  and  most  physical  components  can  be  obtained 
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through  knowledgeable  use  of  the  post-processor  once  the 
calculations  are  complete. 

C.  USA  BOUNDARY  ELEMENT  METHOD  CODE 

The  Underwater  Shock  Analyzer  (USA)  computer  code  [Ref.  6] 
is  a  boundary  element  computer  code  based  on  the  Doubly 
Asymptotic  Approximation  (DAA)  theory  developed  by  Geers  in 
1971  (Ref.  7)  .  Through  the  use  of  the  DAA  theory  and  the 
boundary  element  formulation,  USA  computes  the  acoustic 
pressure  loading  and  added  mass  matrices  which  represent  the 
fluid  surrounding  the  submerged  shell.  The  acoustic  pressure 
loading  and  added  mass  are  applied  at  selected  wetted  nodes. 
This  formulation  has  the  benefit  of  significantly  decreasing 
the  number  of  elements  required  to  model  the  submerged  system 
since  external  water  elements  need  not  be  included  in  the 
calculations.  The  reduced  number  of  elements  requires 
substantially  less  time  and  storage  space  to  obtain  a 
solution . 

Ho’wever,  it  must  be  noted  that  this  code  has  limitations 
v/hich  result  directly  from  the  fundamental  assumptions 
associated  with  the  DAA  theory  (Ref.  6) .  First,  DAA  is  not 
theoretically  appropriate  for  concave  or  multiple  structures 
or  near  surface  problems  involving  convex  bodies.  However, 
studies  show  that  only  results  in  highly  shadowed,  closely 
spaced  areas  or  regions  of  strong  concavity  are  affected. 
Secondly,  DAA  requires  that  the  source  of  the  incident  wave  be 
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sufficiently  removed  from  the  structure  since  it  can  only 
account  for  acoustic  waves  and  not  hydrodynamic  flow. 
Finally,  the  DAA  theory  is  based  on  an  early  time  {high 
frequency)  approximation  coupled  with  a  late  time  (low 
frequency)  approximation.  Therefore,  although  the  DAA 
solution  will  be  very  good  at  early  times  when  the  high 
frequency  approximation  is  dominant  and  at  late  times  when  the 
low  frequency  approximation  is  dominant,  it  can  vary 
significantly  from  the  analytic  or  exact  solution  during 
intermediate  times  when  neither  the  high  or  low  frequency 
solution  is  dominant.  A  detailed  description  of  DAA  theory  is 
provided  by  reference  7 . 

D.  COMPUTER  CODE  VERIFICATION 

Since  the  USA/DYNA3D  interface  was  new  and  had  not  been 
tested,  some  effort  was  expended  on  performing  a  verification 
of  the  performance  of  the  code.  To  perform  the  verification, 
two  cases  with  known  analytic  results  were  modelled  using  the 
USA/DYNA3D  code.  The  first  case  was  a  quarter  cylinder  and 
the  second  was  an  infinite  cylinder.  Results  were 
satisfactory  for  both  cases  and  the  code  interface  is  believed 
to  be  performing  correctly. 

1.  DETAILED  DESCRIPTION  OF  THE  SPHERICAL  MODEL 

The  numerical  study  was  performed  on  a  quarter 
symmetry  model  of  a  sphere  containing  150  elements.  Figure 
II. 1  shows  the  model  and  figure  II. 2  shows  the  test  geometry. 


Figure  II. 1.  Elastic  sphere  test  case  model. 
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Figure  II. 2.  Elastic  cylinder  test  case  geometry. 


The  thickness  to  diameter  ratio  of  the  shell  is  1  to  50  and 
the  shell  is  constructed  of  steel.  The  excitation  is  provided 
by  a  very  small  step  pressure  wave.  As  a  result,  the  shell 
response  is  considered  to  be  completely  elastic.  The  case  was 
run  using  the  elastic  material  model  of  DYNA3D  and,  since 
results  are  being  compared  to  the  analytic  results  found  in 
reference  8,  the  same  material  and  water  properties  as  those 
found  in  reference  8  were  used.  As  stated  in  reference  8,  the 
exact  solution  is  obtained  from  separation  of  variables  as 
shown  in  reference  9.  The  material  and  water  properties  used 
are  listed  below; 

Steel  Properties 

Young's  Modulus  E=206.84  GPa 

Poisson's  ratio  v=0.33 

Mass  density  p=7784.5  kg/m 

Water  Properties 

Sound  speed  c=1461.2  m/s 

Density  p=999.6  kg/m 

The  numerical  results  using  the  US.^/DYN.A3D  combination 
for  the  above  test  case  compare  favorably  with  the  exact 
results.  The  normalized  results  are  shown  in  figure  II. 3.  It 
can  be  seen  that  the  numerical  results  lag  the  exact  results, 
but  the  difference  is  negligible.  Copies  of  the  INGRID  pre¬ 
processor  input  and  output  as  well  as  the  USA  pre-processor 
inputs  are  provided  in  Appendix  A. 
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Figure  11.3.  Spherical  model  verification  results 


2.  DETAILED  DESCRIPTION  OF  THE  INFINITE  CYLINDER  MODEL 

The  infinite  cylinder  model  was  run  using  the  same 
material  and  water  properties  shown  above.  Figure  II. 4  shows 
the  geometry  used  for  the  analyses  and,  as  shown  in  figure 
II. 5,  a  single  ring  of  elements  was  used  to  model  the  infinite 
cylinder  by  enforcing  symmetry  boundary  conditions  on  each  end 
of  the  model.  In  addition,  since  this  is  a  two  dimensional 
problem,  the  TWODIM  option  in  USA  was  used  to  generate  the 
added  mass  and  DAA  matrices.  Further,  the  value  of  the  Ti 
variable  was  set  to  0.0.  r)  is  the  factor  that  accounts  for 
curvature.  This  resulted  in  a  DAAl  solution  for  comparison  to 
a  known  analytic  D.iA.l  solution.  The  first  model  attempted  had 
a  longitudinal  length  of  0.1  inches.  Hov;ever,  it  was 
discovered  that  this  resulted  in  a  oscillatory  solution  as 
shovm  in  the  first  graph  in  Figure  II. 6.  A  similar 
oscillation  occurred  on  the  reverse  side  of  the  cylinder  as 
shown  in  Figure  1 1. 7.  After  a  check  of  the  input  data  to 
ensure  that  the  problem  was  not  caused  by  numerical 
instability,  it  was  hypothesized  that  the  oscillation  was 
caused  by  residual  three  dimensional  effects  caused  by  the 
finite  width  of  the  model.  As  a  test  of  this  hypothesis,  two 
additional  models  were  run  with  widths  of  0.01  and  0.001 
inches.  As  shown  in  Figure  II. 6  and  II.  7,  reduction  in  width 
progressively  reduced  the  oscillations  on  both  the  front  and 
back  of  the  cylinder.  At  0.001  inches,  oscillations  are 
absent . 
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Figure  II. 6,  Infinite  cylinder  oscillation  reduction 
(near  element). 
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Figure  II. 7.  Infinite  cylinder  model  oscillation  removal 
(far  element). 
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The  final  results  from  the  0.001  inch  model  with  a  T| 
variable  value  of  0.0  were  compared  to  the  analytical  exact 
and  analytical  DAAl  solutions  with  favorable  results  as  shown 
in  figure  II. 8.  It  can  be  seen  that  the  results  on  both  the 
front  and  back  sides  of  the  cylinder  lie  very  close  to  the 
analytic  DAAl  solution. 

A  further  investigation  was  conducted  to  determine 
what  value  of  the  T|  variable  would  result  in  the  numerical 
solution  closest  to  the  analytic  modal  solution.  Values  of 
0.0,  0.25,  0.5,  0.75  and  1.0  were  tested.  The  results  are 
shown  in  figures  II. 8  through  11.12.  Review  of  the  results 
show  that  the  value  of  the  T|  variable  that  provides  the 
results  nearest  the  analytical  modal  solution  varies  depending 
on  time  and  position  on  the  cylinder. 

For  the  front  of  the  cylinder,  a  DAA2  variable  value 
of  0.0  gives  results  fairly  close  to  the  analytical  solution 
for  all  times.  However,  values  of  0.25  and  0.5,  although  they 
do  not  match  closely  at  early  times,  match  more  closely  at 
late  times. 

On  the  reverse  side,  a  value  of  0.0  provides  a  result 
very  near  the  analytical  DAAl  solution,  but  varies 
substantially  from  the  analytical  modal  solution.  Values  of 
0.5  and  0.75  provide  results  near  the  analytical  modal 
solution  v;ith  0.75  being  the  best  result. 

Assuming  that  interest  lies  in  late  time  results  over 
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Figure  II. 8.  Infinite  cylinder  results  with  ri=0.0. 
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Fisrure  II.9.  infinite  cylinder  results  with  r|-0.25. 
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Figure  11.10.  Infinite  cylinder  results  with  ri=0.50, 
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Figure  II.  11.  Infinite  cylinder  results  with  r)s0.75. 
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Figure  II.12.  Infinite  cylinder  results  with  r|=1.0. 
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the  entire  cylinder,  the  results  show  that  the  best  overall 
value  of  r|  for  an  infinite  cylinder  lies  between  0.5  and  0.75. 
More  compact  bodies  will  have  best  results  v;ith  higher  T] 
values . 

The  INGRID  pre-processor  input  and  output  for  the  infinite 
cylinder  verification,  as  well  as  the  USA  code  pre-processor 
input  are  provided  in  Appendix  B. 
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III.  EXPERIMENTAL  METHOD  AMD  DESCRIPTION  OF  MODELS 
A.  DESCRIPTION  OF  PHYSICAL  MODEL  DSBD  FOR  THE  FAR  FIELD  STUDY 

The  physical  model  was  an  unstiffened  right  circular 
cylinder  with  the  following  characteristics. 

Dimensions : 

Length  42  inches  (1.067  m) 

Diameter  12  inches  (0.305  m) 

Weight  60.5  pounds  (27.5  Kg) 

Materials : 

Shell  1/4  inch  thick  6061-T6  Aluminum  (0.64  cm) 

End  Plates  1  inch  thick  6061-T6  Aluminum  (2.54  cm) 

The  cylinders  used  for  this  test  'were  constructed  from 
commercially  available  material.  Fabrication  was  performed  at 
the  Naval  Postgraduate  School.  The  end  plates  were  welded  to 
the  shell  using  a  Tungsten  Inert  Gas  (TIG)  process. 

The  6061-T6  aluminum  was  selected  on  the  basis  of  its  high 
strength  and  strain  rate  insensitivity.  The  material 
properties  of  the  aluminum  used  for  the  shell  were  verified 
using  the  MTS  machine  at  the  Naval  Postgraduate  School. 
Results  of  tensile  testing  determined  that  the  material 
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properties  were  near  nominal  with  a  Young's  modulus  of  10800 
ksi  (75.6  GPa)  and  yield  strength  of  43  ksi  (300  MPa). 

B.  UNDERWATER  EXPLOSION  TEST 

The  underwater  explosion  test  was  performed  at  the  Dynamic 
Testing  Incorporated  (DTI)  facilities  in  Rustburg,  Virginia. 
The  facility  is  in  a  quarry  and  the  depth  of  the  water  is 
approximately  130  feet  (39.6  m)  at  the  location  of  the  test. 
As  a  result,  bottom  reflection  was  not  a  factor  in  the  test. 

The  charge  used  for  the  test  was  60  pounds  (27.3  Kg)  of 
HBX-1.  The  peak  pressure  generated  by  the  charge  v;a3  2360 
psig  (16.3  MPa)  v.'hich  was  substantially  lower  than  the 
calculated  peak  pressure  of  2680  psig  (18.5  MPa)  for  the  60 
pound  (27.3  Kg)  charge  at  a  25  foot  (7.62  m)  standoff 
distance.  The  test  charge  was  activated  by  a  radio  control 
device . 

The  test  depth  for  both  the  charge  and  the  cylinder  was  12 
feet  (3.66  m) .  This  depth  allowed  the  bubble  generated  by  the 
explosion  to  vent  at  the  surface  prior  to  encountering  the 
cylinder  and  eliminated  the  possibility  of  a  bubble  pulse.  In 
addition,  the  12  foot  (3.66  m)  depth  provided  a  clear  pressure 
cutof  f  . 

The  cylinder  was  held  in  place  with  a  crane  rig  and  the 
charge  was  suspended  from  a  float.  Distance  and  alignment  of 
the  charge  to  the  cylinder  was  established  and  maintained 
using  a  tensioned  span  wire  from  the  charge  float  to  the 
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cylinder  support  rig.  Post-shot  calculations  found  the 
arrival  time  of  the  shock  wave  to  be  consistent  with  a 
distance  of  25  feet  (3.66  m)  and  sound  of  speed  in  water  of 
4800  ft/sec  (1463  m/s).  Test  profile  and  arrangement  are 
provided  as  figures  III.l  and  III. 2. 

Strain  measurement  was  performed  using  CEA-06-250UW-350 
strain  gages.  These  are  general  purpose  strain  gages  with  an 
optimum  range  of  ±  1500  microstrain  and  are  good  for  both 
static  and  dynamic  test  measurements.  The  strain  gages  were 
bonded  to  the  cylinder  using  a  M  bond  200  by  a  instrumentation 
technician  employed  by  DTI.  All  pre-shot  calibration  and 
connection  were  performed  by  DTI  technicians. 

The  test  called  for  14  total  strain  gages  (seven  to 
measure  hoop  strains  and  seven  to  measure  axial  strains)  .  Of 
the  fourteen  strain  gages,  three  failed.  The  dynamic  range  of 
the  test  exceeded  the  optimum  range  of  the  strain  gages  by  a 
significant  factor.  This  is  the  most  probable  cause  of  the 
high  strain  gage  failure  rate.  The  instrumentation  diagram 
for  the  test  is  provided  as  figure  III. 3.  The  strain  gage 
located  at  B1  was  placed  nearest  the  charge  during  the  test. 
Strain  gage  output  was  filtered  at  2000  Hz.  Locations  noted 
on  figure  III. 3  will  be  used  for  reference  throughout  the 
remainder  of  the  report. 

Slight  damage  to  the  cylinder  was  noted  upon  completion  of 
the  test.  Post-shot  investigations  found  all  strain  gages 
firmly  attached  to  the  cylinder  at  the  locations  specified  in 
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Figure  Iii.l.  undex  test  profile. 
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Figure  III. 2 
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Undex  test  general  arrangement . 
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Figure  III. 3.  Undex  test  instrumentation  diagreun. 
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the  instrumentation  diagram.  However,  some  water  intrusion 
was  noted  under  the  protective  coating  of  several  of  the 
strain  gages.  This  intrusion  may  also  have  played  a  part  in 
the  strain  gage  failures.  The  results  of  the  test  were 
forwarded  to  the  Naval  Postgraduate  School.  A  copy  of  the 
report  is  enclosed  as  Appendix  C. 

C .  NUMERICAL  MODEL 

This  study  was  performed  using  two  primary  mesh  densities. 
The  low  density,  full  model  mesh  (figure  III.  4)  was  used  for 
rotation,  shell  type  and  quadrature  sensitivity  analyses.  The 
high  density  quarter  model  was  used  to  perform  direct 
comparison  to  experimental  results  and  examine  end  effects. 


Figure  III. 4.  Low  density,  full  model 

The  computational  efficiency  of  the  quarter  model  allowed  a 
more  refined  mesh  without  a  subsequent  increase  in 
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computational  time  or  random  access  memory  storage  capability. 
A  sample  quarter  model  was  run  and  results  checked  against  a 
full  model  with  the  same  mesh  configuration  to  certify  that 
the  symmetry  boundary  conditions  used  to  form  the  quarter 
model  were  valid.  The  refined  mesh  quarter  model  is  shown  in 
figure  III. 5. 


Figure  III. 5.  Refined  mesh  model. 

In  addition  to  the  two  model'  noted  above,  several 
additional  quarter  models  with  varying  mesh  density  were  run 
to  verify  mesh  size  independence  of  the  quarter  model  results. 
It  was  found  that  the  most  critical  locations  for  the  mesh 
sensitivity  check  were  the  locations  with  the  highest  strain. 
The  areas  with  the  highest  strain  were  located  near  each  end 
on  the  side  of  the  cylinder  located  nearest  the  explosive 
charge.  Figure  III.  6  shows  the  strain  pattern  on  the  surface 


of  cylinder  side  nearest  the  charge.  The  high  strain 
locations  are  symmetrically  located  16.5  inches  (0.42  m)  from 
the  axial  midpoint  of  the  cylinder.  The  other  region  of 
significant  plastic  strain  was  located  on  the  surface  of  the 
reverse  side  of  the  cylinder  at  the  axial  midpoint.  Figure 
III. 7  shows  the  effective  plastic  strain  pattern  for  this 
location.  Effective  plastic  strain  is  defined  by  the 
relation : 

C  =  f(C  “C  ~e  )^+(e  ~e 

P  3  '•  *'2p'  '  2p  '•3P'  ^^3P  ^ip'  J 


where 

®ip<  ®2p  represent 

the 

principal  plastic 

strain 

components  [Ref.  10]. 

The 

near 

side  high  strain 

regions 

cover  a  much  smaller 

area 

than 

the  reverse  side 

region.  That  is,  much  higher  strain  gradients  occurred  on  the 
near  side  compared  to  other  locations  on  the  cylinder.  This 
condition  plays  a  significant  roll  in  mesh  design  and 
integration  time  increment  selection. 

Figures  III. 8  through  III. 10  show  the  results  of  the  mesh 
sensitivity  test.  It  was  found  that  strains  in  the  axial 
direction  v/ere  more  sensitive  to  mesh  density  than  hoop  strain 
results.  Figure  III. 8  shows  the  strain  at  the  surface  of  the 
cylinder  at  the  point  nearest  the  charge  (location  B1 ) .  This 
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Figure  III. 6.  Effective  plastic  strain  pattern  on  cylinder 
side  nearest  the  explosive  charge. 
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Figure  III. 7.  Effective  plastic  strain  pattern  on  cylinder 
side  most  remote  from  the  charge. 
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location  has  no  permanent  plastic  strain.  It  can  be  seen  that 
there  is  no  significant  difference  between  the  results  for  the 
three  mesh  densities  checked.  Figure  III. 10  shows  strain 
results  for  the  surface  of  the  shell  at  the  point  most  remote 
from  the  charge  in  the  circumferential  direction  at  the  axial 
midplane  (location  B3 ) .  This  location  had  the  second  highest 
strain  of  the  positions  checked.  Although  there  is  a  slight 
difference  between  the  three  different  mesh  results,  it  is 
apparent  that  these  differences  are  insignificant  when 
compared  to  the  overall  plastic  strain.  Figure  III.  9  shows 
the  strain  results  for  the  locations  that  experienced  the 
highest  plastic  strains  (locations  A1  and  Cl) .  The  difference 
in  the  hoop  direction  is  noticeable  but  small  enough  to  be 
neglected.  However,  the  results  in  the  axial  direction  are 
significant  with  a  30  percent  variance  between  the  average 
plastic  strains  for  the  high  density  mesh  and  medium  density 
mesh.  .A.dditional  refinement  was  not  possible  due  to  random 
access  iviemory  limitations  on  the  system,  used  to  perform  the 
analysis.  On  the  basis  of  the  above  results  it  was  determined 
that  the  medium  mesh  model  was  adequate  for  comparison  of 
numerical  to  experimental  results  for  all  hoop  strains  and  all 
axial  strains  except  at  the  locations  near  the  end  on  the  side 
nearest  the  charge.  The  high  density  mesh  was  used  for  the 
axial  strain  comparison  at  the  remaining  locations.  Care  was 
taken  to  ensure  that  the  mesh  was  as  uniform  as  possible  for 
both  the  full  and  the  quarter  mesh  models  to  avoid  problems 
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Figure  III. 8.  Mesh  sensitivity  comparison  for  surface  of 
shell  located  nearest  the  explosive  charge 
(location  Bl) . 
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Figure  III.  10.  Mesh  sensitivity  comparison  for  point  on 
cylinder  most  circiimerentially  remote  from 
the  charge  (location  B3) . 


with  mesh  reflection  as  noted  earlier  in  this  report. 

Thin  shell  elements  were  used  for  both  the  shell  and  end 
plates.  Since  relatively  small  out  of  plane  displacements 
were  encountered  in  the  test  model,  it  was  determined  that  the 
four  node  Belytschko/Lin/Tsay  shell  formulation  [Ref  11] , 
which  is  the  default  formulation  for  VEC/DYNA3D,  was  adequate 
for  the  analysis.  A  Hughes/Liu  [Ref.  12]  shell  model  and  a 
eight  node  brick  shell  model  were  also  run  for  comparison. 

The  Belytschko/Lin/Tsay  shell  was  selected  over  the 
Hughes/Liu  shell  and  8  node  brick  shell  formulation  because  of 
its  higher  relative  computational  efficiency. 

The  aluminum  v/as  treated  as  a  kinematic/isotropic 
elastic/plastic  material  with  no  strain  rate  sensitivity. 
Research  has  shown  that  shock  velocities  much  higher  than  the 
velocities  encountered  in  the  test  are  required  to  induce 
strain  rate  sensitivity  in  6061-T6  aluminum. 

The  pressure  input  for  the  model  was  obtained  from  the 
free  field  pressure  transducer  time  record  of  the  underwater 
explosion  test.  The  17000  point  trace  was  numerically 
condensed  to  100  points  and  entered  into  the  TIMINT  pre¬ 
processor  of  USA  using  the  VARLIN  (variable  linear)  option. 
Figure  III. 11  shows  the  pressure  profile  used  for  the 
analysis.  Free  surface  effects  were  neglected  and  the  speed 
of  sound  in  water  used  for  the  test  was  4800  ft/sec  (1463  m/s) 
since  the  test  was  performed  in  fresh  water  at  approximately 
40  degrees  Fahrenheit  (4.5  degrees  centigrade). 
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IV.  FAR  FIELD  STUDY  RESULTS 


A.  EXPERIMENTAL  TO  NUMERICAL  COMPARISON 

As  described  earlier  in  the  report,  an  underwater 
explosion  test  was  conducted  at  the  Dynamic  Testing, 
Incorporated  facility  in  Rustburg,  Virginia.  The  test 
included  a  side-on  attack  of  a  cylinder  with  a  stand  off 
distance  of  25  feet  (7.62  m)  using  a  60  pound  (27.3  Kg)  HBX-1 
charge.  Fourteen  strain  gages  were  attached  to  the  cylinder, 
of  which  eleven  provided  useable  data.  Four  statements  can  be 
made  about  the  results.  First,  the  numerical  results  compared 
v;ell  with  the  experimental  results  qualitatively.  That  is,  the 
numerical  response  had  the  same  general  shape  as  the 
experimental  results  and  it  predicted  compression  and  tension 
correctly.  There  v;a3  one  exception  to  the  above  statement  at 
position  B3  (Figure  IV. 9).  The  numerical  model  indicated  a 
tensile  a.xial  strain  at  position  B3  v;hile  the  experimental 
data  indicated  a  compressive  strain.  Physically,  it  can  be 
observed  that  the  shock  wave  is  spherical  and  initially 
strikes  the  cylinder  center.  This  places  the  cylinder  in 
bending.  Therefore,  tensile  strain  is  expected  in  the  axial 
direction  on  the  reverse  side  of  the  cylinder.  It  is  believed 
that  the  poles  on  the  axial  strain  gage  at  position  B3  were 
reversed  resulting  in  an  error  in  sign  of  the  data  returned  by 
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numerical  and  experimental  results  match  more  closely  in  areas 
v;ith  lov;er  values  of  total  strain.  Finally,  axial  strains 
were  affected  more  than  hoop  strains.  Charge  size  factors 
were  eliminated  as  a  possible  cause  of  the  magnitude 
differential  since  the  measured  pressure  profile  was  used  to 
perform  the  post  underwater  explosion  test  numerical 
calculations.  In  addition,  the  possibility  of  the  charge  being 
located  closer  to  the  shell  than  the  specified  standoff 
distance  was  eliminated  by  comparing  the  actual  shock  wave 
travel  time  measured  from  the  strain  gage  traces  to  the 
expected  shock  wave  travel  time  calculated  for  the  speed  of 
sound  in  water  for  fresh  water  at  40  degrees  Fahrenheit  (4.4 
degree  centigrade) .  The  results  indicated  less  than  two 
inches  difference  between  the  calculated  and  measured  values 


for  stand  off  distance. 


Third,  the  frequency  of  oscillation  of  the  numerical  data 
v;as  lower  than  the  experimental  data.  The  higher  frequency 
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oscillation  in  the  physical  model  compared  to  the  numerical 
model  indicates  that  the  experimental  model  is  "stiffer"  than 
the  numerical  model.  This  is  an  unexpected  result,  since 
numerical  finite  element  solutions  are  normally  expected  to  be 
stiffer  than  the  physical  model.  In  addition,  the  numerical 
results  for  axial  strain  tended  to  “ring"  at  all  locations. 
The  "ringing"  is  not  a  significant  factor  for  hoop  strains. 
It  should  also  be  noted  that  the  "ringing"  is  heaviest  at  the 
front  and  back  of  the  cylinder  at  the  center.  The  causes  of 
the  "ringing"  and  the  high  stiffness  of  the  physical  model 
have  not  been  determined  and  are  a  topic  of  additional  study. 

Finally,  there  is  an  unexpected  asymmetry  in  the 
experimental  results.  The  axial  strain  gage  at  position  Cl 
{figure  IV. 10)  measured  50%  lower  than  the  axial  strain  gage 
at  Al  (figure  IV. 1)  and  the  hoop  strain  gage  at  position  C2 
(figure  IV. 11)  measured  nearly  50  percent  higher  than  the  hoop 
strain  gage  at  position  A2  (figure  IV. 2).  Failure  of  strain 
gages  at  positions  Al,  Cl,  and  C2  prevented  additional 
comparisons.  The  asymmetric  results  can  result  from  two 
factors.  The  shell  may  have  been  rotated  from  the  expected 
orientation  by  underwater  currents  or  by  forces  placed  on  the 
cylinder  and  rigging  by  the  instrumentation  cables  or  there 
could  have  been  a  failure  in  the  bonding  between  the  strain 
gage  and  the  cylinder  surface  on  one  or  more  strain  gages. 

Figures  IV.  1  through  IV. 11  provide  the  results  of  the 
numerical  to  experimental  data  comparison. 
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Figure  IV. 1.  Experimental /numerical  comparison  for 

position  Al  axial  strain. 


< 

(X 


05 

O 

tr 

o 


TIME  (SEC) 

Figure  IV. 2.  Experimental/numerical  comparison  for 

position  A2  hoop  strain. 
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Experimental/numerical  comparison  for 
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Experimental/numerical  comparison  for 
position  B1  axial  strain. 
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Experimental/numerical  comparison  for 
position  B2  hoop  strain. 


Figure  IV. 6. 
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position  B2  axial  strain. 


TIME  (SEC) 

Experimental /numerical  cougar ison  for 
position  B3  hoop  strain. 


Figure  IV, 8. 


0.5x10'’ 


z 

<  0 

cr 

t— 

w 

o 

tr 

O  A 

f  -0.5x10'’ 


-I.OxlO** 

0  .002  .004  .006  .008  .010 

TIME  (SEC) 

Figure  IV. 11.  Experimental /numerical  comparison  for 
position  C2  hoop  strain. 


B.  SENSITIVITY  ANALYSES 

A  series  of  sensitivity  analyses  were  performed  in  an 
effort  to  explain  the  differences  between  the  numerical  and 
experimental  results  noted  in  the  previous  section.  In 
addition,  these  analyses  provided  additional  insight  into  the 
relative  importance  of  various  factors  in  the  performance  of 
underwater  explosion  tests  and  the  associated  calculations. 
Seven  sensitivity  analyses  were  performed.  The  first  was  the 
mesh  sensitivity  test.  The  results  of  this  analysis  have 
already  been  discussed.  The  other  six  analyses  were,  end 
effect,  shell  element  formulation,  integration  time  increment 
length,  quadrature,  rotational  position  and  T)  value 
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sensitivity  checks.  The  results  of  these  analyses  are 

provided  in  the  following  subsections. 

1.  END  EFFECT  SENSITIVITY  ANALYSIS 

As  previously  noted,  the  most  severe  deformation 
occurred  at  locations  near  the  end  of  the  cylinder  (positions 
A1  and  Cl) .  Two  processes  cause  this  phenomena.  First,  the 
relatively  large  mass  of  the  end  plates  apply  large  inertial 
forces  to  the  cylinder  shell  near  the  end  plates.  Second,  the 
one  inch  thick  end  plates  are  very  stiff  and  their  lack  of 
flexibility  causes  the  weaker  material  of  the  shell  near  the 
end  plates  to  deform  in  response  to  applied  forces.  A 
examination  of  the  numerical  and  experimental  data  reveals 
that  these  effects  are  concentrated  near  the  end  plates  and 
result  in  large  strain  gradients.  This  means  that  elements  on 
either  side  of  a  selected  element  near  the  end  of  the  cylinder 
can  have  significantly  different  strain  values.  Accurate 
placement  of  strain  gages  within  this  region  and  careful  mesh 
design  along  'with  adequately  short  time  integration  increments 
are  critical  in  obtaining  satisfactory  results  in  a  numerical 
to  experimental  data  comparison.  In  addition,  as  stated 
earlier,  the  end  plates  are  attached  to  the  shell  using  a 
Tungsten  Inert  Gas  process.  This  welding  process  results  in 
high  temperatures  near  the  end  of  the  cylinder.  Since  the 
aluminum  for  this  model  is  at  a  peak  hardened  condition,  this 
process  could  result  in  a  change  of  the  material  properties 
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near  the  end  cf  the  cylinder  that  can  only  be  restored  by 
performing  the  age  hardening  process  again  after  the  welding 
is  complete.  These  factors  can  result  in  an  uncertainty  in 
the  expected  strain  compared  to  what  might  occur  under  ideal 
circumstances . 

The  mesh  sensitivity  results  clearly  display  the 
importance  of  mesh  design  within  this  region.  However,  even 
with  proper  design,  the  large  gradients  can  result  in 
significant  differences  between  the  predicted  and  actual 
strains  since  the  strain  computed  for  the  element  is  an 
average  of  the  strain  over  the  entire  element  vice  a  strain  at 
a  specific  point.  The  best  possible  results  would  be  obtained 
in  these  regions  with  large  gradients  if  the  mesh  could  be 
refined  such  that  the  size  of  the  elements  is  the  same  size  as 
the  gage  length  of  the  strain  gage.  However,  this  would 
result  in  a  prohibitively  large  number  of  elements  and  a 
subsequent  increase  in  problem  solution  times.  These  problems 
can  be  overcome  by  placing  strain  gages  in  areas  that  are 
expected  to  have  consistently  increasing  or  decreasing  strains 
and  then  ensuring  that  the  mesh  is  designed  so  that  the  strain 
gage  location  is  at  the  center  of  the  element.  If  possible, 
large  gradient  regions  should  be  avoided.  If  strain  gages 
must  be  placed  in  a  high  gradient  region,  then  the  strain 
gages  should  be  placed  to  one  side  or  the  other  of  the  minimum 
or  maximum  strain  location.  Placement  at  the  minimum  or 
maximum  point  will  result  in  an  error  since  the  average  for 
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the  element  v;ill  lie  above  a  minimum  or  below  a  maximum  if  the 
element  is  not  the  same  size  as  the  gage  length  of  the  strain 
gage. 

In  this  study,  the  strain  gages  located  at  Al  and  Cl 
were  located  at  the  point  of  highest  compressive  strain. 
Therefore,  a  study  was  performed  to  determine  the  relative 
importance  of  the  noted  location  factors.  Figures  IV.  12 
through  IV. 16  show  the  results  of  this  study.  Strains  of  two 
additional  elements  nearer  the  end  were  compared  to  the 
measured  strain  and  the  actual  strain  gage  location.  Elements 
nearer  the  end  plate  were  selected  since  the  welding  effects 


Figure  IV. 12.  End  effect  sensitivity  results. 
(Al  axial) 
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Figure  IV. 13.  End  effect  sensitivity  results. 
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Figure  IV.  14.  End  effect  sensitivity  results. 

(A2  Axial) 
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Figure  IV. 15.  End  effect  sensitivity  results 
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Figure  IV. 16.  End  effect  sensitivity  results 

(C2  Hoop) 


described  in  the  previous  paragraph  would  tend  to  move  the 
high  strain  location  nearer  to  the  end  plate  by  weakening  the 
material  near  the  end  plate.  Only  the  positions  with  useable 
experimental  results  are  shown.  In  four  of  the  five  cases 
(positions  Al  axial,  A2  hoop.  Cl  axial  and  C2  hoop),  if 
asymmetry  effects  are  taken  into  account,  the  element  one 
nearer  to  the  end  from  the  actual  strain  location  provides  a 
better  estimate  of  the  actual  strain  measured  during  the 
underwater  explosion  test.  At  the  fifth  location  (A2  axial), 
the  second  element  closer  to  the  end  provides  the  best 
results.  These  results  require  additional  study  to  separate 
and  quantify  the  effect  of  the  phenomena. 

2.  SHELL  FORMULATION,  QUADRATURE  RULE  AND  INTEGRATION 
TIME  INCREMENT  SENSITIVITY  ANALYSES 
In  addition  to  the  above  end  effects,  there  was  some 
concern  that  the  mid  plane  reference  for  the  thin  shell 
element  would  result  in  a  greater  flexible  length  than  the 
actual  physical  model.  This  concern  was  based  on  the  fact 
that  the  mass  and  stiffness  of  the  end  plates  is  concentrated 
into  a  planar  surface  co-located  with  the  mid  plane  of  the  end 
plate  in  the  thin  shell  analysis.  This  resulted  in  the  shell 
portion  of  the  structure  being  one  inch  longer  in  the 
numerical  model  than  the  physical  model.  This  problem  could 
have  been  avoided  by  using  the  Hughes/Liu  foruvulation  and 
shifting  the  reference  plane  to  the  inner  surface  of  the 


shell.  To  resolve  this  issue  a  study  was  conducted  to  compare 
the  performance  of  different  types  of  thin  and  thick  shell 
element  formulations. 

Results  from  the  Belytschko/Tsay/Lin  shell  formulation 
were  compared  to  results  from  the  same  model  using  the 
Hughes/Liu  shell  formulation.  As  stated  earlier,  the 
Belytschko/Lin/Tsay  shell  has  the  advantages  of  increased 
computational  efficiency  and  a  high  degree  of  stability  with 
large  deformations  at  the  expense  of  reduced  accuracy  at  high 
levels  of  plastic  strain.  The  major  difference  betv;een  the 
two  formulation  stems  from  the  fact  that  the  element  normal 
direction  is  updated  periodically  in  the  Hughes/Liu 
formulation.  The  Belytschko/Lin/Tsay  formulation  assumes 
negligible  cut  of  plane  deformations,  and  therefore,  does  not 
update  the  shell  normal.  As  a  result,  the  inaccuracy  of  the 
Belytschko/Lin/Tsay  formulation  will  increase  as  shell 
deformation  becomes  significant. 

The  models  used  to  compare  the  two  formulations  were 
identical  in  all  aspects  with  the  exception  of  the  shell 
formulation.  The  center  line  plane  was  used  for  the  reference 
on  both  models.  The  results  confirmed  that  the  strain  levels 
encountered  in  this  underwater  explosion  test  were  small 
enough  to  support  use  of  the  Belytschko/Lin/Tsay  formulation. 
However,  it  was  apparent  that  differences  did  occur  for 
positions  with  significant  plastic  strain  in  the  axial 
direction  (Positions  A1 ,  A2,  B3 ,  Cl,  and  C2 ) .  Although  the 
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differences  in  these  cases  were  not  significant  enough  to 
require  use  of  the  Hughes/Liu  formulation,  it  is  also  noted 
that  higher  strain  may  result  in  larger  differences. 
Therefore  the  Belytschko/Lin/Tsay  formulation  should  not  be 
used  in  cases  where  significant  denting  occurs  unless 
stability  problems  occur  while  using  the  Hughes/Liu 
formulation . 

As  stated  earlier,  the  presence  of  high  strain 
gradients  near  the  end  plates  causes  small  changes  in  end 
condition  or  distance  to  be  significant.  When  it  became 
apparent  that  end  effects  would  be  important  in  the  results  an 
investigation  was  performed  to  determine  if  an  eight  node 
brick  shell  formulation  would  provide  more  accurate  results 
near  the  end  of  the  cylinder.  The  thin  shell  formulation 
results  as  well  as  the  experimental  results  were  compared  to 
results  from  a  model  computed  using  eight  node  brick  shell 
elements.  All  three  formulations  are  compared  to  experimental 
results  in  figures  IV. 17  through  IV. 27.  The  following 
information  can  be  gleaned  from  the  plotted  results.  First, 
it  is  apparent  that  the  greatest  differences  occur  near  the 
positions  with  the  highest  strains.  At  the  same  time,  it  can 
be  noted  that  there  is  virtually  no  difference  at  the 
locations  with  no  permanent  strain.  Second,  as  shown  in 
figures  IV. 17,  IV. 18,  IV. 24,  IV. 25  and  IV. 27,  it  is  clear  that 
there  is  a  significant  difference  between  the  eight  node  brick 
shell  results  and  the  Belytschko/Lin/Tsay  results  at  the 
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locations  with  high  levels  of  permanent  plastic  strain. 
However,  contrary  to  the  expected  results,  the  eight  node 
brick  shell  results  move  further  from  the  expected  values  than 
the  other  formulations.  It  is  also  noted  that  the  Hughes/Liu 
formulation  lies  between  the  eight  node  brick  shell  and  the 
Belytschko/Lin/Tsay  formulation. 

Additional  research  was  performed  to  determine  the 
cause  of  the  disparities.  The  study  revealed  that  the  eight 
node  brick  shell  is  sensitive  to  integration  time  increment 
and  v/ill  move  marginally  closer  to  the  thin  shell  results  if 
time  integration  is  cut  in  half.  However,  the  overall  shift 
is  only  about  10  percent  of  the  total  difference.  Quadrature 
rule  (number  of  points  used  in  the  Gauss  quadrature  numerical 
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Figure  IV. 17.  Shell  formulation  sensitivity  results. 
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Figure  IV. 18.  Shell  formulation  sensitivity  results. 
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Figure  IV. 19.  Shell  formulation  sensitivity  results 
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ire  IV. 20.  Shell  formulation  sensitivity  results. 
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Figure  IV. 21.  Shell  formulation  sensitivity  results. 
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ire  IV. 22.  Shell  formulation  sensitivity  results. 
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Figure  IV. 23.  Shell  formulation  sensitivity  results. 
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ire  IV. 24. 
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Figure  IV. 25.  Shell  formulation  sensitivity  results. 
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ire  IV. 26.  Shell  formulation  sensitivity  results. 
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Figure  IV. 27.  Shell  formulation  sensitivity  results. 
(C2  Hoop) 
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integration  scheme)  proved  to  be  a  more  significant  effect. 
Use  of  five  point  quadrature  moved  the  three  results  closer 
together  while  having  the  most  profound  effect  on  the 
Belytschko/Lin/Tsay  formulation.  Again,  the  affect  of 
quadrature  rule  affected  the  thick  shell  results  only 
marginally.  Figure  IV. 28  shows  the  combined  results  for  the 
location  with  greatest  plastic  deformation. 


Figure  IV. 28.  Effect  of  changing  quadrature  rule  and  tine 
integration  increment  at  location  of  highest 
strain  (Al  and  Cl  hoop) . 

In  summary,  all  three  formulations  appear  to  be 
satisfactory  as  long  as  care  is  used  in  designing  the  mesh  and 
selecting  the  integration  time  and  quadrature  integration 
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rule  parameters .  Specif ically,  when  using  Belytschko/Lin/Tsay 
formulations  in  areas  with  relatively  high  strain,  the  number 
of  quadrature  points  should  be  increased  until  stable  results 
are  achieved.  When  using  eight  node  brick  shell  elements, 
integration  time  increment  must  be  selected  with  care  but 
number  of  quadrature  points  seems  to  be  less  critical.  The 
Hughes/Liu  formulation  appeared  to  be  relatively  insensitive 
to  both  quadrature  rule  and  integration  time  increment. 

Reference  13  provides  some  useful  thumbrules  for 
selection  of  time  increments.  The  following  criteria  are 
recommended . 

At  »  0.9- -  for  brick  shells 

At  »  0.9-^  for  thin  shells 
Dc 

V  -  element  volume 
At  -  time  increment 
Ag  -  maximum  surface  area 
D  -  maximum  diagonal 
c  -  speed  of  sound  in  the  material 
Aj,  -  maximum  area  of  any  surface 


The  above  criteria  were  found  to  be  adequate  except  for  the 
highest  strain  areas  where  the  thick  shell  element  rule  did 
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not  provide  stable  solutions.  In  areas  such  as  A1  and  Cl,  a 
value  of  the  integration  time  increment  half  of  the  above 
recommendation  proved  to  be  satisfactory  for  the  eight  node 
brick  shell. 

3.  ROTATION  SENSITIVITY  ANALYSIS 

A  sensitivity  analysis  was  performed  to  determine  the 
effect  of  an  in  plane  rotation  away  from  the  expected 
symmetric  orientation  in  an  effort  to  explain  the  cause  of  the 
asymmetric  results  of  the  underwater  explosion  test.  It  was 
hypothesized  that  an  unplanned  rotation  greater  than  ten 
degrees  v;ould  have  been  detected  by  the  personnel  performing 
the  test.  Four  different  models  were  run  within  this  range 
representing  rotations  of  0.0,  2.5,  5.0  and  10.0  degrees.  The 
results  are  shov.-n  along  with  e.xperimental  results  vjhere 
available  in  figures  IV. 29  through  IV. 42.  The  following 
observations  are  made  concerning  the  results.  First,  the  most 
dramatic  affects  are  on  the  reverse  side  of  the  cylinder  at 
position  B3  (Figures  IV. 37  and  IV. 38).  The  results  show  that 
the  differential  between  the  numerical  and  experimental 
results  at  position  B3  can  be  explained  by  a  six  to  eight 
degree  rotation  from  the  symmetric  configuration.  Rotational 
effects  at  locations  Bl  and  B2  on  the  centerline  (Figures 
IV. 33  through  IV. 36)  are  insignificant.  Hoop  strain  at 
position  C2  (Figure  IV. 41)  is  approximately  60  percent 
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ire  IV. 29.  Rotation  sensitivity  results. 
(Al  Hoop) 
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Figure  IV. 30.  Rotation  sensitivity  results. 

(Al  Axial) 
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ire  IV. 31.  Rotation  sensitivity  results. 
(A2  Hoop) 
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Figure  IV. 32.  Rotation  sensitivity  results. 
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ire  IV. 35.  Rotation  sensitivity  results. 
(B2  Hoop) 
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Figure  IV. 36.  Rotation  sensitivity  results. 

(B2  Axial) 
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jire  IV. 37.  Rotation  sensitivity  results 
(B3  Hoop) 
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rigure  IV. 38.  Rotation  sensitivity  results 
(B3  Axial) 
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Figure  IV. 39.  Rotation  sensitivity  results. 
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Figure  IV. 40.  Rotation  sensitivity  results. 
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ire  IV. 41.  Rotation  sensitivity  results 
(C2  Hoop) 
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Figure  IV. 42.  Rotation  sensitivity  results 
(C2  Axial) 


higher  than  the  hoop  strain  at  position  A2  (Figure  IV. 32) 
with  a  rotation  of  ten  degrees.  This  is  also  consistent  with 
the  experimental  data.  Similar  positive  results  were  obtained 
for  positions  A1  and  Cl  axial  strains.  It  was  .further 
discovered  that  rotating  the  cylinder  about  its  axis  could 
further  improve  the  results.  However,  even  though  these 
rotations  did  improve  the  results,  significant  differences 
still  exist  between  the  experimental  and  numerical  strains  at 
the  ends  of  the  cylinders.  Although  it  is  clear  that  the 
model  can  account  for  rotational  effects,  it  is  also  clear 
that  other  factors  are  causing  the  large  differences.  Once 
again,  welding  affects  are  suspected  to  be  the  probable  cause. 

The  important  point  to  note  out  of  these  results  is 
that  even  small  rotations  from  expected  orientation  can  result 
in  significant  errors  on  in  expected  results.  Therefore 
extreme  care  must  be  taken  to  ensure  that  instrumentation 
cable  tension  and  other  unanticipated  factors  do  not  cause 
undetected  rotations. 

4.  PHYSICAL  FINDINGS 

a.  RESPONSE  MODES 

It  was  determined  that  a  cylinder  subject  to  a  side 
on  explosion  will  have  three  primary  response  modes.  The 
first  mode  is  an  accordion  motion.  The  accordion  motion 
results  from  the  compression  and  subsequent  release  of  the 
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cylinder  in  the  axial  direction.  Figure  IV. 43  shows  a  plot  of 
points  located  at  the  center  of  each  end  plate.  It  is  clear 


TIME  (sec) 

Figure  IV. 43.  Cylinder  accordion  motion. 

that  the  two  end  plates  are  travelling  in  opposite  directions 
at  the  same  time  generating  the  accordion  motion. 

The  cylinder  is  also  subject  to  a  whipping  mode 
parallel  to  the  direction  of  shock  wave  travel.  The  whipping 
mode  is  the  most  significant  motion  experienced  by  the 
cylinder  and  is  caused  as  a  result  of  the  curvature  of  the 
shock  wave.  In  the  symmetric  situation,  the  shock  v/ave  will 
come  in  contact  with  the  center  of  the  cylinder  first.  This 
will  cause  the  center  to  move  first,  followed  by  the  ends. 
The  cylinder  will  then  move  in  an  oscillatory  motion  that  is 
a  function  of  the  stiffness  and  mass  distribution  of  the 
cylinder  and  the  'water  surrounding  the  cylinder.  Fj.gure  IV. 44 
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Figure  IV. 44.  Cylinder  whipping  motion  in  plane  parallel 
to  shock  wave  direction. 

shov/s  a  plot  of  a  points  located  at  the  center  and  ends  of  a 
line  located  parallel  to  the  axis  on  the  near  side  of  the 
cylinder.  The  plot  shows  that  the  end  plates  are  moving  in 
the  opposite  direction  of  the  cylinder  throughout  the 
transient  response  of  the  cylinder.  Figure  17.45  shows  a 
scale  factor  20  drav.’ing  of  the  cylinder  at  two  different 
times.  The  cylinder's  opposite  direction  of  curvature  at  the 
two  different  times  is  a  result  of  the  whipping  motion. 

The  final  response  mode  noted  was  a  breathing 
motion  in  the  plane  perpendicular  to  the  shock  wave  direction 
of  travel.  Although  breathing  motion  also  occurred  in  the 
direction  parallel  to  the  shock  wave  travel,  it  was  not  as 
obvious  since  the  miuch  larger  v;hipping  motion  turned  out  to  be 
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TIME  ( sec ) 

Figure  IV. 46.  Cylinder  breathing  motion  perpendicular  to 
the  shock  wave  direction  of  travel. 

point  is  moving  in  a  direction  opposite  to  the  lower  point 
throughout  the  transient  response  of  the  cylinder.  The 
breathing  motion  is  also  caused  by  the  compression  and 
subsequent  release  of  the  cylinder.  Figure  IV. 47  provides  an 
illustration  of  the  breathing  motion.  The  two  scale  factor  40 
drav;ings  are  for  two  separate  times  and  show  the  shell  first 
bowed  inward  to'ward  the  axis  and  then  outward  away  from  the 
axis . 


b.  ROTATIONAL  EFFECTS 

Plastic  strain  fringe  plots  generated  as  a  result 
of  the  rotation  sensitivity  analysis  revealed  some  interesting 
information  on  the  causes  of  the  strain  distribution 
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Figure  IV. 47.  Illustration  of  cylinder  breathing  mode  at 
two  different  times  (scale  factor  40) . 

experienced  by  the  cylinder.  The  experimental  results 
included  a  reduction  in  the  strain  at  the  rear  of  the  cylinder 
at  position  B3 ,  a  decrease  at  A2  relative  to  C2  and  an 
increase  at  Al  relative  to  Cl.  The  fringe  plots  show  why  this 


strain  distribution  occurs.  Figure  IV. 48  shows  the  effective 
plastic  strain  distribution  for  a  7.5  degree  rotation.  The 
left  side  of  the  cylinder  is  nearest  the  charge.  The  results 
show  that  the  rotation  tends  to  diffuse  the  strain  around  the 
cylinder  on  the  near  end  while  concentrating  it  at  the  far 
end.  This  causes  the  distribution  noted  for  positions  Al,  Cl, 
A2  and  C2 .  At  the  same  time,  the  high  stress  region  on  the 
reverse  side  of  the  cylinder  tends  to  move  away  from  the 
charge.  This  placed  location  B3  in  a  lower  strain  region 
which  led  to  the  experimental  and  numerical  results  noted  at 
position  B3 . 
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Figure  IV. 48.  Effective  plastic  strain  distribution  for 
near  and  remote  side  of  cylinder  with  7.5 
degrees  of  rotation. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS. 

1 .  NUMERICAL  MODELLING 

Two  general  conclusions  can  be  reached  from  the 
material  contained  in  this  report.  First,  the  USA/DYNA3d 
connection  is  successful  and  can  replicate  the  response  of 
simple  analytical  models. 

Second,  numerical  modeling  can  predict  the  response  of 
a  simple  cylinder  to  an  underwater  explosion.  Far  field 
numerical  predictions  generally  match  experimental  results  if 
rotation  and  end  effects  resulting  from  fabrication  caused 
material  property  changes  are  correctly  modelled.  It  was 
found  that  results  in  high  strain  areas  are  extremely 
sensitive  to  shell  formulation,  mesh  design,  quadrature  rule 
and  integration  time  increment.  The  best  results  were  achieve 
with  brick  s.hell  elements.  However,  the  eight  node  brick 
shell  required  substantially  longer  computation  times  to 
achieve  the  desired  results  because  of  the  need  to  reduce  the 
integration  time  increment.  In  addition,  it  was  found  that 
thin  shell  formulations  can  also  provide  correct  results. 
However,  results  for  the  Belytschko/Lin/Tsay  formulation 
appear  to  be  very  sensitive  to  the  number  of  quadrature  points 
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used  for  the  numerical  integration  scheme  in  high  strain 
areas . 

2.  PHYSICAL  ASPECTS 

Primary  damage  areas  are  near  the  ends  of  the  cylinder 
on  the  side  nearest  the  charge  where  the  stiff,  heavy,  flat 
end  plates  caused  a  concentration  of  the  effective  plastic 
strain.  Damage  also  occurred  on  the  reverse  side  as  a  result 
of  a  bending  effect  similar  to  that  described  in  the  near 
field  results.  The  cylinder  experienced  breathing,  whipping 
and  accordion  response  modes. 

In  addition,  it  was  discovered  chat  rotation  tends  to 
diffuse  strain  on  the  end  nearest  the  charge  while 
concentrating  the  strain  at  the  far  end  on  the  side  nearest 
the  charge.  The  high  strain  area  located  at  the  center  of  the 
cylinder  on  the  reverse  side  tends  to  migrate  toward  the  end 
most  remote  from  the  charge. 

B .  RECOMMENDATIONS . 

1.  TOPICS  FOR  ADDITIONAL  STUDY. 

a.  WELDING  FABRICATION  EFFECTS. 

An  analysis  should  be  performed  to  quantify  the 
relative  effect  that  the  change  in  material  properties 


generated 

by  the 

welding  fabrication 

process 

has 

on  the 

numerical 

results 

This  analysis 

could 

include  the 

measurement  of  material  properties  near 

a  weldment . 

These 

81 


properties  could  then  be  added  as  a  separate  material  in  the 
numerical  model . 

b.  EIGHT  NODE  BRICK  SHELL  SENSITIVITY  ANALYSES. 

Although  it  was  fairly  clear  that  the  eight  node 
brick  shell  formulation  comes  closest  to  predicting  the 
overall  response  of  the  shell,  it  was  also  noted  that  the 
formulation  is  very  sensitive  to  integration  time  increment  in 
areas  with  high  strain.  Commonly  used  thumbrules  did  not 
appear  to  be  adequate  in  this  case.  In  addition,  additional 
analyses  need  to  be  performed  to  determine  the  mesh 
sensitivity  of  the  eight  node  brick  shell  in  this  model. 

c.  FAILURE  CRITERIA. 

This  study  was  performed  on  a  model  with  relatively 
lov;  total  plastic  strain  (less  than  one  percent)  .  In  order 
to  deal  with  larger  strains,  a  failure  model  must  be 
introduced  into  the  material  modelling  of  the  cylinder.  The 
model  should  include  structural  instability  as  well  as 
material  rupture  criteria. 

d.  NEAR  FIELD  EXPERIMENTATION 

Although  the  numerical  predictions  appear  to  be 
physically  correct,  the  physical  results  obtained  using  them 
cannot  be  assumed  completely  correct  until  they  are  confirmed 
with  experimental  results.  A  study  should  be  conducted  to 
compare  near  field  experimental  results  with  numerical 
predictions . 
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2.  RECOMMENDATIONS  TO  IMPROVE  TEST  CONTROL. 

Several  factors  made  the  comparison  of  the  numerical 
to  experimental  results  difficult.  If  properly  controlled, 
the  analysis  process  could  be  simplified.  First,  rotation  of 
the  cylinder  must  be  carefully  controlled.  Second,  unless 
specifically  required,  high  strain  gradient  areas  should  be 
avoided.  Placement  of  the  strain  gages  becomes  critical  in 
these  locations  as  does  mesh  design  and  integration  time 
increment.  If  these  areas  cannot  be  avoided,  additional 
sensitivity  analyses  may  be  required  to  determine  the  adequacy 
of  the  mesh  and  integration  time  increment.  Finally,  analysis 
near  welded  seams  should  be  avoided  unless  the  effects  can  be 
quantified.  If  near  weld  analysis  cannot  be  avoided, 
consideration  should  be  given  to  restoring  the  heat  treatment 
after  the  v;eld  process. 
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A.  INGRID  PRE-PROCESSOR  INPUT  DATA 
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B.  INGRID  PRE-PROCESSOR  ODTPDT  DATA 
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c. 


FLUMAS  PRE-PROCESSOR  DATA 


FLUMAS  DATA  FOR  SPHERICAL  SHELL  QUARTER  MODEL 


sph. flu 

sph . geo 

sph.nom  sph.daa  $ 

FLUNAM 

GEONAM 

GRDNAM 

DAANAM 

T  T  F 

T 

$ 

PRTGMT 

PRTTRN 

PRTAMF 

CALCAM 

T  F  F 

T 

$ 

EIGMAF 

TWODIM 

HAFMOD 

QUAMOD 

F  F  T 

T 

$ 

PCHCDS 

NAS TAM 

STOMAS 

STOINV 

F  F  F 

F 

$ 

FRWTFL 

FRWTGE 

FRWTGR 

FRESUR 

F  T  F 

F 

$ 

RENUMB 

STOGMT 

ROTGEO 

ROTQUA 

F  F  T 

F 

$ 

PRTCOE 

STRMAS 

SPHERE 

ROTSYM 

F  F  F 

F 

$ 

OCTMOD 

CAVFLU 

FRWTFV 

EFUDGE 

DYNA 

$ 

MAINKY 

0  171 

0 

150 

$ 

NSTRC 

NSTRF 

NGEN 

NGENF 

0  0  0 

$ 

NBRA 

NCYL 

NCAV 

1.  1. 

$ 

RHO 

CEE 

10 

$ 

NVEC 

1.  1. 

1  . 

1 . 

$ 

CQ(1) 

CQ{2) 

CQ(3) 

CQ(4) 

1 

$ 

NSRADI 

1 

0 

1 


1.  1  150  1 


$  RADI  RAD2 
$  NSORDR 
$  SPHRAD 


JBEG 


JEND  JINC 
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D.  AUGMAT  PRE-PROCESSOR  DATA 

AUGMAT  DATA  FOR  SPHERICAL  SHELL  QUARTER  MODEL 


sph.nom  sph.flu 

sph.geo  sph.pre  $ 

STRNAM 

FLUNAM 

GEONAM 

PRENAM 

F  F  F  F 

$ 

FRWTGE 

FRWTST 

FRWTFL 

PLNWAV 

F  F  T  T 

$ 

FLUSKY 

DAAFRM 

SYMCON 

DOFTAB 

F  F  F  T 

$ 

PRTGMT 

PRTTRN 

PRTSTF 

PRTAUG 

F  F  F 

$ 

MODTRN 

STRLCL 

INTWAT 

DYNA 

$ 

MAINKY 

1  . 

$ 

DAA2 

171  513  3 

3 

$ 

NSTR 

NSFR 

NFRE 

NFTR 

1 

$ 

NSETLC 

0  1  150  1 

$ 

NDICOS 

JSTART 

JSTOP 

JINC 

4 

$ 

NUMCON 

1  1  121 

6 

$ 

ICON 

NSTART 

NSTOP 

NINC 

2  16 

1 

$ 

ICON 

NSTART 

NSTOP 

NINC 

2  37  42 

1 

$ 

ICON 

NSTART 

NSTOP 

NINC 

2  127  167 

5 

$ 

ICON 

NSTART 

NSTOP 

NINC 
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E.  TIMINT  PRE -PROCESSOR  INPUT  DATA 


TIMINT  DATA  FOR  SPHERICAL  SHELL  QUARTER  MODEL 


sph.pre  sph.pos 

$ 

PRENAM 

POSNAM 

sph . rst 

$ 

RESNAM 

0.  0.  10000. 

$ 

XC 

YC 

ZC 

0.  0.  1. 

$ 

sx 

SY 

SZ 

F  F  F  F 

$ 

EXPWAV 

SPLINE 

VARLIN 

PACKET 

F  F  F  F 

$ 

HYPERB 

BUBPUL 

REFSEC 

EXPLOS 

2 

$ 

JPHIST 

1.  0. 

$ 

PNORM 

50. 

$ 

DTHIST 

l.E-6  l.E-6 

$ 

PHIST 

1  0 

$ 

NT  I  NT 

NCHGAL 

0.  .00025 

$ 

STRTIM 

DELTIM 

1000  1000 

$ 

NSAVER 

NRESET 

NSODFL 

0  0  0  0 

$ 

LOCBEG 

LOCRES 

LOCWRT 

NSTART 

F 

$ 

DISPLA 
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A.  INGRID  PRE-PROCESSOR  INPOT  DATA 

INFINITE  CYLINDER  MODEL 

dn3d  vec  term  5.0  plti  0.01  prti  1000. 

mat  1  type  1  e  9.813e+l  pr  0.3  ro  7.85 
shell  quad  0  thick  0.010  endmat 

led  1  2  0.0  0.0  10.0  0.0 

led  2  2  0.0  0.000025  10.0  0.000025 


plan  3 


0  0  0 

0 

-1  0 

0.00001 

symm 

0  0  -.0005 

0  0 

-1  0.001 

symm 

0  0  .0005 

0  0 

1  0.001 

symm 

start 

-1  6 

-11 

/ 

-1  6 

-11 

t 

1  2  ; 
-1  .  0  . 

1 . 

-1.  0. 

1 . 

-.0005 

.0005 

alio 

3  3 

0 

3  1.0 

d  0  1  0 

0  2 

0 

pri  -1 

-3  ; 

-1 

-3  ; 

;  1  -1.0 

O 

o 

mate  1 
end 

end 
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B.  INGRID  PRE-PROCESSOR  OUTPUT  DATA 
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C.  FLUMAS  PRE- PROCESSOR  INPUT  DATA 


FLUMAS  DATA  FOR 
cyl.flu  cyl.geo 
T  T  T  T 
T  T  T  F 
F  F  T  T 
F  F  F  F 
F  T  F  F 
F  F  F  F 
F  F  F  F 
DYNA 

0  42  0  20 
0  0  0 

1.  1. 


20  1 


INFINITE  CYLINDER  MODEL 

cyl.nom  cyl.daa  $  FLUNAM  GEONAM  GRDNAM  DAANAM 

$  PRTGMT  PRTTRN  PRTAMF  CALCAM 

$  EIGMAF  TWODIM  HAFMOD  QUAMOD 

$  PCHCDS  NASTAM  STOMAS  STOINV 

$  FRWTFL  FRWTGE  FRWTGR  FRESUR 

$  RENUMB  STOGMT  ROTGEO  ROTQUA 

$  PRTCOE  STRMAS  SPHERE  ROTSYM 

$  OCTMOD  CAVFLU  FRWTFV  EFUDGE 

$  MAINKY 

$  NSTRC  NSTRF  NGEN  NGENF 

$  NBRA  NCYL  NCAV 

$  RHO  CEE 

$  NVEC 

$  CQ(1)  CQ ( 2 ) 

$  DHALF  CXHF  CYHF  CZHF 

$  NSRADI 

$  RADI  RAD2  JBEG  JEND  JI 
$  NSORDR 

$  NORD  JBEG  JEND  JINC 


NGEN 

NCAV 


JBEG  JEND  JINC 
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D.  AUGMAT  PR2- PROCESSOR  DATA 


AUGMAT  DATA  FOR 

INFINITE 

CYLINDER 

MODEL 

cyl .nom  cyl .flu 

cyl . geo 

cyl.pre  $ 

STRNAM 

FLUNAM 

GEONAM 

F  F  F  F 

$ 

FRWTGE 

FRWTST 

FRWTFL 

F  F  T  T 

$ 

FLUSKY 

DAAFRM 

SYMCON 

T  T  F  T 

$ 

PRTGMT 

PRTTRN 

PRTSTF 

F  F  F 

$ 

MOCTRN 

STRLCL 

INTWAT 

DYNA 

$ 

MAINKY 

0.75 

$ 

DAA2 

42  126  3  3 

$ 

NSTR 

NSFR 

NFRE 

1 

$ 

NSETLC 

0  1  20  1 

$ 

NDICOS 

JSTART 

JSTOP 

2 

$ 

NUMCON 

2  17  6 

$ 

ICON 

NSTART 

NSTOP 

2  33  38  5 

$ 

ICON 

NSTART 

NSTOP 

PRENAM 

PLNWAV 

DOFTAB 

PRTAUG 


NFTR 

JINC 

NINC 

NINC 
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B.  TIMINT  PRE-PROCESSOR  DATA 


TIME  INTEGRATION 

DATA  FOR 

INFINITE  CYLINDER  MODEL 

cyl .pre 

cyl .pos 

$ 

PRENAM 

POSNAM 

cyl . rst 

$ 

RESNAM 

10000  . 

0.0000 

0.000000 

$ 

XC 

YC 

ZC 

1.0000 

0.0000 

0 . 000000 

$ 

SX 

SY 

SZ 

F  F  F  F 

$ 

EXPWAV 

SPLINE 

VARLIN 

PACKET 

F  F  F  F 

$ 

HYPERB 

BUB PUL 

REFSEC 

EXPLOS 

2 

$ 

J PHI ST 

1.  0. 

$ 

PNORM 

HYDPRE 

50  . 

$ 

DTHIST 

l.e-6  1. 

e-6 

$ 

PHISTd)  PHIST{2) 

1  0 

$ 

NTINT 

NCHGAL 

0.  0.00025 

$ 

STRTIM 

DELTIM 

10000  10000 

$ 

NSAVER 

NRESET 

NSODFL 

0  0  0  0 

$ 

LOCBEG 

LOCRES 

LOCWRT 

NS TART 

F 

$ 

DISPLA 
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DYNAMIC  TESTING,  INC. 


P.O.  Box  494  0  Rustburg,  Virginia  24SBB-0494  •  (B04)  B46  0244  •  Fax  (B04)  B46-2197 


DTI  R  TSG 
DTI91-183 
4  September  1991 


United  States  Naval  Postgraduate  School 
Department  of  Mechanical  Engineering 
Monterey.  Calitornia  93943-5100 

Attention:  Dr.  Young  Shin 

Subject:  Cylinder  Test  Results,  DTI  Job  No.  131 

Enclosure:  (1)  Instrumentation  Time- Histories,  Test  No.  1 

(2)  Instrumentation  Time-Histories,  Test  No.  2 

(3)  Test  Photographs 

(4)  Data  Provided  in  ASCII  Format 

Gentlemen: 

Dynamic  Testing.  Inc.  (DTI|.  is  pleased  to  provide  the  foliowing  results  for  the  cylinder 
tests  conducted  at  our  facility  on  20  and  22  August  1991. 

The  first  test,  conducted  on  Tuesday,  20  August  1991,  using  a  60-pound  HBX 
cylindrical  charge,  consisted  of  one  cylinder  oriented  for  an  end-on  shot  with  a  28-toot 
standoff  at  a  depth  of  12  feet.  There  were  a  total  of  7  strain  gauges  per  axis  for  a 
total  of  14  strain  gauges.  One  free-field  pressure  gauge  was  located  28  feet  from  the 
charge,  but  was  positioned  in  such  a  way  that  no  reflection  from  the  cylinder  occurred. 
There  were  three  pair  of  strain  gauges  located  in  the  center  of  the  cylinder,  one  at  Bl, 
B2,  and  B3.  Two  pair  of  strain  gauges  were  located  4.5  inches  from  each  end 
oriented  at  Al,  A3,  Cl,  and  C3.  Before  the  test  It  was  noted  that  gauge  AlC  was 
not  working  properly,  so  we  oriented  the  cylinder  so  that  the  end,  designated  C,  was 
toward  the  blast  (see  Figure  1). 

The  strain  gauges  were  oriented  to  monitor  longitudinal  and  circumferential  strains. 
Predicted  analysis  Indicated  a  .2  percent  strain  at  30  feet.  No  visual  damage  was 
observed. 

The  second  test,  conducted  Thursday,  22  August  1991,  using  a  60-pound  HBX 
cylindrical  charge,  consisted  of  one  cylinder  oriented  for  a  side  shot  with  a  25-foot 
standoff  at  a  depth  of  12  feet.  There  were  a  total  of  7  strain  gauges  per  axis  for  a 
total  of  14  strain  gauges.  One  free-field  pressure  gauge  was  located  25  feet  from  the 
charge,  and  positioned  so  no  reflection  from  the  cylinder  occurred.  There  were  three 
pair  of  strain  gauges  located  in  the  center  of  the  cylinder,  one  at  Bl,  B2,  and  B3. 
Two  pair  of  strain  gauges  were  located  4.5  inches  from  each  end  oriented  at  Al.  A2. 
Cl.  and  C2  (see  above). 

The  strain  gauges  were  oriented  to  monitor  longitudinal  and  circumferential  strains. 
Predicted  analysis  indicated  a  .48  percent  strain  at  25  feet.  No  visual  damage  was 
observed. 


Dr.  Young  Sltin 

Uiiiied  Stales  Naval  Postgraduate  School 
DTI91  183 
4  September  1991 
Page  2 


The  test  configuration  consisted  of  one  cylinder  being  suspended  by  the  crane  and  one 
chaigc  suspended  from  a  float  (see  Figures  2a  and  2b  for  test  coidiguialions). 

We  trust  you  will  find  the  enclosed  satisfactory.  It  was  a  pleasure  woiking  with  you 
on  this  project. 

If  you  have  any  questions  or  comments,  please  do  not  hesitate  to  contact  Mr.  Randy 
Faittield  or  the  undeisigned. 

Sincerely. 

DYNAMIC  TESTING.  INC. 

Tony  S.  Grigsby 
Inslrumenlatiou  Technician 

TSG:mdl 

cc:  G.G.  Amir 

G.D.  Snyder 
R.D.  Faittield 
W.G.  Lyon 


DTIOl-183 


Dr.  Young  Shin 

United  States  Naval  Postgraduate  School 
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Figure  2a.  Test  Configuration  for  Test  No.  1 
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Figure  2b.  Test  Configuration  for  Test  No.  2 
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